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ABSTRACT* 

Investigations have been performed aimed at developing plastic nuclear track 

detectors into precise and quantitative tools for recording and measuring 

multicharged, heavy particles. Accurate track etch rate measurements as a 

function of LET have been performed for cellulose nitrate and Lexan plastic 

detectors. This was done using a variety of incident charged particle types 

and energies. The effect of aging of latent tracks. in Lexan in different 

gaseous atmospheres has been investigated. Range distributions of high 
14 

energy N particle Bevatron beams in nuclear emulsion were measured. 
Investigation of charge resolution and Bragg peak measurements were carried 
out using plastic nuclear track detectors. 


*The track geometry portion of this report was previously published in 
Nuclear Instruments and Methods 97_ (1971) 483-489. 



INTRODUCTION 


The purpose of this effort was to carry out investigations aimed at 

1 2 

refining the techniques and methods of plastic nuclear track detectors. 9 
The goal is to develop plastic nuclear track detectors into precise and 
quantitative tools for heavy particle measurements. Some of the areas 
investigated include: The behavior of the track etch rate as a function 

of particle LET in Lexan and cellulose nitrate, the response of different 
types of cellulose nitrate as well as of the different surfaces of a 
single sheet of the same material, the response as a function of the etching 
time, and the aging effect of the latent tracks as a function of storage 
atmosphere. An analysis of the response of cellulose nitrate and Lexan 
is also given. Finally, the detectors were utilized in practical applications 
involving the test of charge resolution and the measurement of the Bragg 
peak of multi charged, high energy particle beams from the Bevatron. 

CELLULOSE NITRATE 

A stack of cellulose nitrate layers of the type USF #4 was exposed to a 
1 6 

stopping beam of 0 particles at the Berkeley Bevatron. The individual 
detector layers were etched at 40°C in 6.25 N NaOH solution without stirring. 
Two etch times were used, 5.0 and 20.0 hours. Measurements of the track 
etch rate, Vy, were performed on both sets of layers and for both surfaces 
of each detector sheet. The results are shown in Figures 1, 2 and 3. Here 
V-j. in microns per hour is plotted against the particle residual range in 
microns. It is observed that on a log - log plot, the data is well- 
approximated by a straight line in each of the cases. 
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In Figure 1 is shown the Vy measurements for both the "air" and the "metal" 

surface of the detector. This refers to the casting process for preparing 

cellulose nitrate where a nitrocellulose solution is poured into a metal pan 

and is then allowed to dry. It is observed that for a 5.0 hour etch, the data 

from both surfaces join smoothly together. 

The Vj data for the 20.0 hr. etch is given in Figures 2 and 3. Inter- 

comparisons with the data for the 5.0 etch shows that the "air surface" data 

for the 20.0 hr. etch fits well together with both sets of data for the 5.0 

etch. However, the "metal surface" data fo,r the 20.0 hr. etch is measurably 

different from the other three sets. This difference suggests a sensitivity 

difference between the two plastic surfaces. Apparently, there is a gradual 

surface-to-surface transition in the properties of the cast material which 

becomes more pronounced with increasing etch times. These observations were 

further substantiated by an even longer etch time of 100 hr. Here, large 

differences in the track structure between the two surfaces was noted. 

The track etch rate as a function of residual range was also measured 

in USF #3 cellulose nitrate, using samples exposed at the PPA to high energy 
20 

beams of Ne particles. This detector is similar to USF #4 type, however, 
is slightly less sensitive. The track etch rates were measured on one layer. 
The tracks due to the particles were then traced to their stopping points to 
find the residual range that applied for the measured etch rate. A total of 
48 tracks were traced through, and four layers were required to reach the 
stopping points of all of the particles. Because of the non-uniform thickness 
of the CN layers, the thickness of each layer at the point of passage of each 
particle had to be measured. 

The results are shown in Figure 4. The track etch rate, V y, is plotted 

as a function of the particle residual range, R. The etch rate calibration 

obtained by drawing a straight line on this log - log plot through the top 

20 

group of tracks which are assumed to be Ne particles is 
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V T = 3083 LET 3 ' 626 , 

where \lj is in p/hour and LET is LET^q in MeV/p. Using this calibration, the 

dashed curves shown were drawn for the isotopes indicated. 

It can be seen that the measured points fall nicely into two groups. The 

logical charge assignment for the upper and most populated group is the nominal 
20 

beam ion, Ne. From this assumption and the calibration curve drawn through 

this group of points, the dashed curves follow. It is assumed that the ions 

represented by the other groups of points are produced through nuclear inter- 

20 

actions of the incident Ne beam. For this reason it is not surprising to 
find some of the points not falling on the dashed line corresponding to the 
most abundant isotope in nature. 

The range of track etch rates measureable in this case is largely governed 
by the constant etch time of the detector used for these measurements. A 
longer etch time would have allowed the measurement of much lower etch rates, 
and a shorter etch time would have allowed the measurement of somewhat greater 
etch rates. The same calibration would not necessarily follow, however. This 
is because 1) the etch rate calibration can only be expected to be represented 
by a power function in LET over a limited range of LET values, and 2) the etch 
rate is a function of the etch time because it varies with the depth in the 
sample. 


- 3 - 




R (u) 

Figure 1 


- 4 - 


(M/hr) 



R (u) 


Figure 2 
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LEXAN 

It is well known that latent tracks in plastic nuclear track detectors 
a 

are affected by various environmental parameters such as UV light, oxygen in 
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the atmosphere, etc. A study was performed to determine the affect of 

storage of latent tracks in Lexan stored in various gaseous atmospheres 

including 0 2 , N 2 Ar, C0 2 and vacuum. Samples of Lexan exposed to 10.1 MeV/ 

40 

nucleon Ar ions were stored in the various atmospheres, in the dark, and 
then processed at various intervals of time. The results are shown in Figure 
5. As is observed, samples stored in inert atmospheres of N 2 , C0 2 and Ar 
showed no measurable change in \lj. Samples stored in vacuum showed a 
definite decrease in with storage time. A sample stored in air at 0°C 
in the dark showed no measurable increase in even after a 1,000 day 
storage period. * 

1 g 

Measurements of UV-irradiated Lexan were made using stopping 0 
particle beam of the Bevatron. An UV-irradiation of 32 hours on each 
side (16-32-16 hours per side) was carried out using a 1,000 watt GE 
mercury arc lamp at a distance of 12 cm from the arc. The results are 
given in Figures 6, 7 and 8. It is observed that again the data are 
well -approximated by a straight line on a log - log plot. In this case, 
all the data join smoothly together. While a greater variation is 
observed in the response over a single surface, as compared with, for 
example, USF #4 cellulose nitrate, a much smaller variation in response 
is observed from one surface of the detector to the next. It should be 
pointed out that some of the scattering of the data is the result of 
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measurements of tracks due to secondaries stopping in the detector which 
were produced by the interaction of the primary beam with the stopping 
material. 
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COMPARISON OF RESPONSE: CELLULOSE NITRATE AND LEXAN 


Plastic nuclear track detectors such as cellulose nitrate, UV irradiated 
Lexan, and non-irradiated Lexan differ not only in sensitivity but also in the 
manner in which the etch rate, Vy, varies as a function of particle LET. Two 
detectors may be comparable in sensitivity but because of the differences in 
the response would be useful in the measurement of different portions of the 
particle LET spectrum. In this section, the response curves are compared for 
one type of cellulose nitrate, USF#3, and UV irradiated as well as non-irradiated 
Lexan. 

The basic expression of the etch rate response function of a detector is 
the track etch rate as a function of the particle LET 35Q . Here it is assumed 
that 350 eV is the proper value w for Vy(Z,A,E) = Vy(LET^) where LET u (Z,A,E). 

A somewhat better parameter than the track etch rate is the normalized track 
etch rate, 


v = 


for cellulose nitrate and non-UV irradiated Lexan or 

V, 

v = 


'1 


V 


G1 


for UV irradiated Lexan. Here V and V rn are the etch rates corrected for UV 

1 G.1 

attenuation. The quantity v is more useful than Vy because the ratio of Vy/V^ 
is the governing factor for both for particle registration (v > esc 6) and 
measurability of Vy in the region v » 1. 

In Figure 9 are shown measurements of v(LET 3 g Q ) for several different plastic 
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processing combinations. The two curves of greatest interest are the USF 
#3 curve and the UV-irradiated Lexan curve. They cover the greatest range 
in v(l < v < 200). Several comparisons can be made between these two response 
curves: 

1. For all values of v, 4 ] n is greater for CN than for 

d ln LET 350 

Lexan. This means the CN gives a much more sensitive measure- 
ment of LET. Thus a given percentage accuracy measurement of 
v gives a more accurate value of LET in CN than in Lexan. 

2. For the measurable range of v (0 < v < 200) the range of 
LET 35 Q values for Lexan (0.1 < LET 35Q < 1.0) is much greater 
than the values for CN. This means that CN is only useful 
over a much more limited range of LET-^q values. In fact, 
as can be seen by the non - UV-irradiated curve, Lexan 
without UV irradiation can be used to measure greater values 
of LET 350 yet. (If the curve is extrapolated, the upper 
limit of measurable v ^200, is reached at LET 35Q - 2.5 
MeV/y). 


3. 


d^ In v ^ 

(d In LET 350 ) 


<0 CN 
>0 Lexan 


This positive curvature of Lexan and negative curvature of CN 
on the log - log plot implies that CN yields its greatest LET 
measurement accuracy at low LET and v values, and Lexan its 
greatest accuracy at high LET values. It should be mentioned, 
however, that the shift in the Lexan curve at about v = 20 is 
probably due to the accumulation of etch products in the track 
becoming important for v £ 20. The effect of this saturation is 
variable, resulting in an actual decrease rather than the 
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apparent increase in LET measurement accuracy. Except for 
this variability in the Vy(LET 350 ) curve, accuracy with Lexan 
for v £ 2 is comparable to that of CN but just covers a 
different range of LET. 

4. The critical values of LET 35Q , LET 35Q (v = 1) = LET 35Q c of 

Lexan and CN are approximately the same. That is, in a sense, 
UV-irradiated Lexan is as sensitive as CN. This is only true, 
however, for <5 = 90° tracks. . At 6 = 90°, a particle will 
register to the same vaiue of LET 35Q in UV-irradiated Lexan 
and CN. The value of LET 33 q cannot be directly compared, 
however, because the stopping properties of the two materials 
are different. In CN the critical value of 3 is still greater 
than in Lexan, indicating a greater sensitivity. 
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COMPARISON OF LEXAN AND CN IN TERMS OF USE 

The following discussion still refers to Figure 9 and the remarks of 
the previous discussion still apply. 

1. Z determination 

Both Lexan and CN are useful for Z determination but each has 
its greatest use in a particular situation. It should be 
remarked that the basis for all of the discussion in this 
paper is the assumption of the consistency of the calibration. 

If tracks of a given LET have very nonreproducible values of 
v, of course the discussion based on the calibration curve is 
not true for the real situation. The parameter R reg is 
maximum residual range at which a particle can register. 

a. Particles for which both Vy and R are known. 

In this case a single measurement of Vy(R) gives Z 

much more accurately in CN (about 5 times more accurately). 

b. V-j-(R) measured over entire non-etched-through portion of R . 

In this situation CN is clearly superior. (1) Fewer values 
of Vy(R) must be measured. (2) Even with fewer values of 
Vy(R) measured, at 6 = 90° CN gives ^ 3.6 times the accuracy 
of Lexan in Z determination. (3) With decreasing dip angles 
the approximate accuracy of Lexan remains the same, but for 
CN it increases as (sin 6)'°'^. (4) The fraction of bulk 

material to be removed to obtain the greatest accuracy in Z 
is 0.50 for Lexan but only 0.27 for CN, allowing greater 
accuracy in < 5 , and a measurements in CN. 

The above discussion assumes equality in the detector thickness for Lexan 
and CN. 
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c. Measurements of Z in a relatively thin stack («R ) where 

- - - .... 3 .... 

the particle stopping points are random . 

In this situation Lexan has a slight advantage over CN in 
one respect. Because of the much greater range of LET values 
measurable, a larger fraction of the registered particles 
will produce measured Z values than in CN. Many of the 
tracks in CN will be etched through holes. However, this 
disadvantage may be offset by the greater stopping density. 
Thus, the number/cm^ which are not etched through may be 
the same in CN as in Lexan. Furthermore, for those tracks 
which are not etched through in CN, the accuracy in Z 
measurement will be better than in Lexan. 

Direct measement of the LET spectrum . 

Direct measurement of the LET spectrum means measuring V T and 
thus LET for each particle track. For this use Lexan is clearly 
superior because of the greater range of LET that is measurable. 
In CN only the range 0.1 < LET 35Q <0.3 could be measured and 
even in this case several samples with different etch times 
would have to be employed. For Lexan the entire range of 0.1 < 
LET 350 < 4.0 could be measured by using different etch times and 
processing with and without UV. For a single etch time with UV 
the range of LET 33 q <1. The only advantages offered by CN are 
better geometry tracks and greater accuracy in the LET measured, 
but the accuracy of Lexan is quite sufficient. 



Measurement of particle fluences, LET spectra by track counting , 

and the “density" parts of the stopping density measurements . 

Here CN is clearly superior for several reasons. 

a. CN is more sensitive (also greater density, p, hence 
greater stopping power per volume) giving greater fluences 
and better counting statistics. 

b. CN is more isotropic in its response. It does not 
discriminate as strongly as Lexan against the low dip, 

6, tracks. Thus when the beam is not isotropic, CN 
would enable the values of the fluences, LET spectra, 
and stopping densities as averaged over solid angle 

to be more closely approximated. 

c. The cutoff in etch rate as R is approached is much 
more rapid in CN. This means that the measured track 
numbers would be much less subject to scanning efficiency 
variations. A much smaller fraction of the total number 
of tracks would be of the marginal character. 
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HEAVY PARTICLE MEASUREMENTS 

In order to test the usefulness of plastic nuclear track detectors in 
practical situations some measurements were performed utilizing the high 
energy, ^250 MeV/nucleon, particle beams now available at the Berkeley 
Bevatron. First, however, nuclear emulsions were utilized in order to 
define the particle beam characteristics. 

14 

1 . Ranges of N Ions in G. 5 Emulsion 

The ranges of 99 tracks were measured in an Ilford G.5 emulsion stack 
exposed at the Bevatron. Several points along each track were measured to 
obtain the "true" range as well as the projected range. All tracks showing 
an interaction were systematically excluded. 

The distribution in projected range, R , for all 99 tracks is shown in 

r 

Figure 10. On this scale the distribution in true range, R, would appear 
essentially the same and has not been plotted. The peak region (3.90 < R,R 

r 

< 4.15 cm) is plotted on a greatly expanded scale for R p in Figure 11 and R 
in Figure 12. 

The information contained in. the range distributions is the central 
value and the moments about this value. For our purposes, the quantities 
of interest are the central value and the RMS spread as given by the standard 
deviation. For gaussian distributions the best estimate of the central value 
is the arithmetic average. If all 99 measured tracks are included we obtain 
the following average values for R and R p 

R = 4.004 ± 0.043 cm 
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R = 3.998 ± 0.043 cm 

r 


and the standard deviations of a single measurements are 


= 0.430 ± 0.031 cm 



=0.430 ± 0.031 cm. 


Examination of the particle range distributions shows, however, that 
the above results suffer because of the several ranges which deviate 
appreciably from the main body of distribution. Only 11 of 99 tracks or 
11.1% rather than the expected 31.7% of the particles have ranges greater 
than one standard deviation from the average. This implies that the 
distributions are not gaussian and the measured dispersion is contributed 
mostly by the few largely deviating tracks and does not accurately represent 
the dispersion of the peak. Thus a better average can be obtained by 
restricting ourselves to the tracks in the immediate vicinity of the peak. 
For the 84 tracks in the range 3.90 < R,R p < 4.15 cm the average ranges are: 

R = 4.0308 ± 0.0039 cm 


Rp = 4.0244 ± 0.0037 cm 

and the standard deviations of a single measurement are 

= 0.036 ± 0.003 cm 

a D = 0.0338 ± 0.0026 cm 
R P 

Careful- examination shows that even these restricted portions of the distributions 
are non-gaussian. Thus, the final approach to obtaining the true central values 
of the peaks and the stopping beam alone unaltered by nuclear reactions is 
to fit a gaussian curve through the counts in the three "bins" nearest the 
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peak. This is accomplished not by taking the average and standard deviation 
but by fitting 


AN, - <f), AR - exp 


1 - 6 . 


j 2 l 


2o 


to the three measured values of AN. by adjusting the total number, N, the 
average value, R, and the standard deviation, a. This approach gives: 


and 


R = 4.0237 ± 0.0018 cm, 
R p = 4.0200 ± 0.0014 cm, 

a R = 0.0117 ± 0.0013, 


a n = 0.0091 ± 0.0010. 


This implies an energy of 242.0 ± 0.1 MeV/amu at the entrance face of the 
emulsion. 

It can be seen that the counts in the three central bins around the 
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peak do represent the stopping point distribution due to the N ions from 
an almost mono-energetic source which have not undergone nuclear interactions 
or any other range dispersing interactions other than the range straggling 
due to the randomness of the randomness of the energy transfers to electrons. 
For protons of 242 MeV/amu the standard deviation in range due to range 
straggling is 1.10%. Since range straggling goes inversely as the square 
root of the atomic number, this implies a 0.295% straggle for The 

values obtained from the three central bins are: 


100 o R 

= 0.29 ± 0.03% 

R 
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and 


TOO cr^ 

2- = 0.23 ± 0.02% 

R P 

The' difference between the average "true" range and projected range is 
also of interest. It can be interpreted in terms of the multiple scattering 
Using the average obtained from all 99 tracks 

R - R = 0.006 ± 0.061 cm. 

P 

The averages from the restricted range interval, 3.90 < R,Rp < 4.15 cm., 
gives: 

R - R = 0.0064 ±0.0054 cm. 

The gaussian fits to the 3 central bin results give 

* - V °- 0037 * °- 0023 

Each of these differences implies that to obtain a statistically meaningful 
difference between the true length and the projected length, many more 
tracks would have to be taken. Taking tracks individually alleviates this 
problem by removing the spread due to the energy dispersion and the range 

straggle. For all tracks with 3.90 < R,R < 4.15) except numbers 67, 17, 

r 

and 13 which had unusually large values of R - Rp) the average value of 
R - R is given by - 

r 

R"^~R = 0.00442 ± 0.00028 cm. 

r 

and the standard deviation of a single measurement by 

a R - Rp = 0.00254 ± 0.00020" 

Since Rp = R cos 0 , where 0 is the angle made by the track with the average 
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track direction. 


R - R p = R (1 - cos 0 ) - R[1 

W = 2 ( — - P ) = 0.00220 ± 0 

R 


The RMS value of this angle is 


-(1- f)] 

.00014. 


v^= 0.0469 ± 0.0015 radians 
= 2.69 ± 0.09 degrees. 



So 
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2 . Charge Measurements Using Cellulose Nitrate . 

Layers of cellulose nitrate, USF #3, detectors were assembled into stacks 

and exposed normally to the incident ^280 MeV/nucleon, ion beam. The 

incident particle beam was stopped using a 28 layer CN stack (with individual 

layer thickness of approximately 270 microns) after having been degraded in 

2 

energy by a 10.6 g/cm stack of Lexan. The Lexan degrader thickness was 
based on the initial expectation of a 250 MeV/nucleon energy beam. Thus 
only a few straggling particles rather than the expected Bragg peak were 
stopped in the CN layers. 

The detector layers were individually etched at 40°C in 100 ml of 6.25N NaOH 
The layer used for locating the stopping point for the particle was etched for 
5 hours. The immediately adjacent layer was etched for 2.5 hours, the next 
for 30 hours, and the last for 48 hours. 

The stack configuration and etched track parameters are shown in Figure 
13, which is not to scale. The average residual range for each track, R, is 
the distance from the stopping point to the center of the etch cone, L. For 
each track the track etch rate, Vy, is given by 



where t is the etch time for that layer, 
e 

For all the tracks measured Vy has been plotted as a function of R, as 

shown in Figure 14. Two groups of points are apparent. One curve has been 

14 

drawn through the measured points attributed to N particles. The second 

TO 

group of points represents C particles (some beam contamination was present). 
A point from an independent experiment utilizing 10 MeV/nucleon 16 0 particles 
is also shown as is the bulk etch rate, V r = B/t . The bulk etch rate repre- 
sents the smallest value of Vy possible. From Figure 14 it is noted that Vy 
is a function of the etch time. This implies a depth dependence of Vy in the 


detector. 
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As can be seen from Figure 14, the measured response of the CN detector 
is such that charge identification of the CNO group of stopping particles is 
achieved. This charge resolution should be compared with the apparently better 
resolution indicated in Figures 2, 3 and 4. In fact, the tightness of the 
points about the Tines in Figures 2 and 3 implies isotropic separation might 
be possible, even for individual particles. The greater dispersion of Figure 
14 combining the data from several layers, etched for different lengths of 
time. This corresponds to the usual situation in which absolute measurements 
of Z are to be made using calibration data obtained from other sheets of 
plastic. The smaller dispersion in Figures 2, 3 and 4 indicates that in 
cases where single sheets of plastic can be internally calibrated, or where 
single sheets are used for differential measurements, greater accuracy can 
be obtained than with mixed groups of layers. 
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As can be seen from Figure 14, the measured response of the CN detector 
is such that charge identification of the CNO group of stopping particles is 
achieved. This charge resolution should be compared with the apparently better 
resolution indicated in Figures 2, 3 and 4. In fact, the tightness of the 
points about the lines in Figures 2 and 3 implies isotropic separation might 
be possible, even for individual particles. The greater dispersion of Figure 
14 combining the data from several layers, etched for different lengths of 
time. This corresponds to the usual situation in which absolute measurements 
of Z are to be made using calibration data obtained from other sheets of 
plastic. The smaller dispersion in Figures 2, 3 and 4 indicates that in 
cases where single sheets of plastic can be internally calibrated, or where 
single sheets are used for differential measurements, greater accuracy can 
be obtained than with mixed groups of layers. 
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Figure 14 




3.. Bragg Peak Measurement . 

A stack of 30 Daicell (red) cellulose nitrate plastics, each 250 microns 
thick was exposed normally to an incident beam of °0 particles. Enough Lexan 
plastic degrader was added in front of the stack so as to stop the beam roughly 
in the center of the CN stack. An exposure consisting of a single pulse of 
31,154 counts was made. The plastics were then etched for 8.0 hr. at 40°C in 
6.25N NaOH solution. 

Individual layers were then optically scanned and all stopping particle 
tracks recorded. The results are shown. in Figure 15. The width of the primary 

peak is apparently primarily the result of the particle range straggling with 

! • • ' 

the initial momentum spread of the particle beam making very little contribution. 

The subsequent stopping particles are apparently lower charge fragments created 
in interactions. The experiment clearly showed that plastic detectors are 
quite useful in mapping the three-dimensional structure of the Bragg peak of 
heavy particle beams. 
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Figure Captions 

1 fi 

Figure 1. Calibration plot of V T vs. R(y) of . o ions in USF 4 cellulose 

nitrate for a 5.0 hr. etch at 40°C in 6.25N NaOH. Data for 

both "air" and "metal" surfaces is given. 

Figure 2. Calibration plot of V T vs. R(y) of ^0 ions in USF 4 cellulose 

nitrate for a 20.0 hr. etch at 40°C in 6.25N NaOH. "Air" surface 

data only. 

Figure 3. Same as in Figure 2 except the data is for the "metal" surface. 
Note the difference from the previous three sets of data. 

— 20 

Figure 4. Calibration plot of \lj vs. R(y) of Ne ions in USF 3 cellulose 
nitrate. Data is from a single surface. 

Figure 5. Effect of aging of latent tracks in Lexan as a function of aging 
atmospheres. 

Figure 6. Calibration plot of Vy vs. R"(y) of ^0 ions in UV-irradiated 
Lexan processed for 5.0 hr. in a saturated solution at 40°C. 

Figure 7. Calibration plot of Vy vs. R^(y) of ^0 ions in UV-irradiated 
Lexan processed for 20.0 hr. in a saturated solution at 40°C. 

Figure 8. Same as Figure 7 except processed for 100.0 hr. 

Figure 9. Etch rate ratio Vy/Vg as a function of LET 350 for USF 3 cellulose 
nitrate, UV-irradiated Lexan, non-UV-irradiated Lexan. 

Figure 10. Distribution of projected ranges of 250 MeV/nucleon ^N particles 
in Ilford G.5 emulsion. 

Figure 11. Same as Figure 10 except the main peak has been enlarged. 

Figure 12. Distribution of the true ranges. 

Figure 13. Schematic drawing of a processed plastic detector stack. 

Figure 14. Charge resolution in USF 3 cellulose nitrate. 

1 

Figure 15. The measured Bragg peak of stopping 0 particles in cellulose 
nitrate. 
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